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Abstract 

Neutron-deuteron elastic scattering cross sections are calculated at dif- 
ferent energies using modern nucleon-nucleon interactions and the Tucson- 
Melbourne three-nucleon force adjusted to the triton binding energy. Predic- 
tions based on NN forces only underestimate nucleon-deuteron data in the 
minima at higher energies starting around 60 MeV. Adding the three-nucleon 

forces fills up those minima and reduces the discrepancies significantly. 
PACS numbers: 21.30.+y, 21.45,+v, 24.10.-i, 25.10,+s 



1 



Typeset using REVTpjX 



Substantial progress has been made recently in the study of the three-nucleon (3N) sys- 
tem both experimentally and theoretically. The set of data is being significantly enriched 
for cross sections and spin observables in elastic neutron-deuteron (nd) and proton-deuteron 
(pd) scattering and in the 3N breakup process. Theoretical formulations and numerical algo- 
rithms have been significantly improved with the result that 3N bound and scattering states 
can be solved exactly. In addition the nucleon-nucleon (NN) system is still very intensively 
investigated and the increased data set provides a sound foundation for reliable modern 
phase-shift analysis [|lj . Based on these phases modern NN forces have been constructed by 
different groups @-§J. These interactions reproduce the NN data set with unprecedented 
accuracy as measured by a x 2 /datum very close to 1. Although those forces are not yet 
linked to the underlying QCD due to well known reasons, they cover a wide spectrum of 
expected properties. There are local, weakly nonlocal and strongly nonlocal versions, purely 
phenomenological ones and semi-phenomenological ones guided by meson theory. They form 
an interesting basis to study few-nucleon systems. Thus theoretical tools and data are avail- 
able to approach the basic questions about the dynamics of three interacting nucleons: is the 
unperturbed free two-nucleon (2N) interaction predominant or is it significantly modified 
due to the presence of an additional nucleon? Such a modification would be one out of many 
other possible three-nucleon force (3NF) mechanisms. What is their significance altogether? 
If the free NN force is sufficient to a very large extent, are specific NN interactions favored 
in the description of the 3N system? In the future QCD should provide theoretically consis- 
tent NN and 3N forces and specifically the relative importance of the latter ones for binding 
energies and scattering matrices. First steps on that ground are being done in chiral per- 
turbation theory Q. Despite of that still restricted theoretical insight from QCD one can 
go ahead and compare the theoretical predictions obtained with modern NN interactions 
and model 3N forces to experimental 3N data. There might be a clearcut signal coming 
from certain observables which deny to be explained by 3N Hamiltonians based on modern 
NN interactions only. Such a "smoking gun" observable would then put limits on present 
day 3N force models and would be also of great importance to test the future QCD based 
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dynamics. 

The three-nucleon binding energy by itself is a first signature. The modern NN inter- 
actions underbind 3 if, but to a different extent ||. The essentially local ones lack binding 
energy of about 800 keV out of 8.48 MeV whereas the nonlocal CD Bonn interaction Q 
underbinds only by ^500 keV. That nonlocality is generated by the Dirac structure of the 
spin 1/2 nucleon. One might conjecture that more dynamical nonlocalities caused by the 
composite nature of the nucleons might have additional effects, which are not sufficiently 
well mimicked by the variety of present day NN force models. Anyhow, that information 
from 3 H on insufficient dynamics based on present day NN forces only should be enriched 
by further evidence from the 3N continuum. 

Such a search for 3N continuum observables, which could serve as a "smoking gun" for 
3NF effects has been persued ever since 3N continuum calculations have become feasible [|7| . 
With the advent of the optimally tuned NN forces and the feasibility to also include 3NF's 
into 3N continuum calculations the conclusive power of such calculations has increased 
tremendously. It is the aim of this article to point to such a "smoking gun" in the 3N 
continuum based on modern 3N Faddeev calculation. 

Before coming to that let us briefly describe the situation in 3N continuum studies. A 
detailed overview has been given recently |§. The bulk of 3N scattering observables below 
about 100 MeV nucleon lab energy can be described quite well in the NN force picture 
only. A beautiful example is the total nd cross section ||. This most simple picture is 
also quite stable in the sense that the most modern phase-equivalent NN force models yield 
essentially the same predictions. But there are exceptions, "time dependent ones", which 
were removed by subsequent measurements ||10|| , but more important true ones, where the 



data are reconfirmed by independent measurements. Such a distinguished case is the low 
energy vector analyzing power A y in elastic Nd scattering |Tl[| . A drastic discrepancy between 
the predictions based on NN forces only and both nd and pd data has been found. Present 
day 3NF models have insignificant effects and do not remove that discrepancy. It is known, 
that A y depends very sensitively on the 3 Pj NN forces. Thus a trivial explanation might 



be that the 3 P NN phase-shift parameters from modern phase-shift analysis have not been 
settled to the true ones ||13|| . Presently it is an unsolved puzzle. If the reason does not lie in 
the NN forces a 3NF of still unknown properties will be responsible. 

Another possible signature for 3NF effects is the space star configuration in the 3N 
breakup process at 13 MeV |]1U|| . In that configuration the three nucleons emerge from the 
interaction region with equal energies and interparticle angles of 120° in the 3N cm. system. 
Also the plane spanned by the three momenta is perpendicular to the beam direction. Two 
nd measurements agree essentially with each other but deviate from theory (in the NN 
picture only and including 3NF models). The situation poses even more questions since pd 
data deviate very severely from the nd data pointing to unexpectedly large Coulomb force 



effects [14 



In the present study we investigate the angular distribution in elastic Nd scattering, the 
most simple elastic Nd scattering observable. The transition amplitude for this process is 
composed of the nucleon exchange part (PGq 1 ), the direct action of a 3NF and a part having 
its origin in the multiple interactions of 3 nucleons through 2N and 3N forces: 

< <f/\U\<f> > = < <P'\PG 1 + Vi 1] (l + P) + PT + Vf } (1 + P)G Q f\<f> > (1) 

That rescattering part is expressed in terms of a T operator which sums up all multiple 
scattering contributions through the integral equation [12[ 



f\(j) > = tP\(j> > +(1 + tG )vi l \l +P)|0 > +tPG f\ ( f ) > 

+ (l+tG )V} 1 \l + P)G o f\<j>> (2) 

Here Go is the free 3N propagator, t the NN t-matrix, and P the sum of a cyclical and 
anticyclical permutation of three objects. The 3NF V4 is split into 3 parts 
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V 4 = $> 4 W (3) 

4 = 1 



where each one is symmetrical under exchange of two particles. For the n — n exchange 
3NF for instance fl5| this corresponds to the three possible choices of the nucleon, which 



undergoes the (off-shell) tc — N scattering. The asymptotic state \<f> > (\(f)' >) is a product 
of the deuteron wave function and the momentum eigenstate of the third particle. 

The exchange part comprises two processes where the incoming nucleon ends up as a 
constituent of the final deuteron and the constituents of the initial deuteron are free in the 
final state. Due to the nature of this term its contribution to the elastic scattering cross 
section is peaked at backward angles. The contribution from the driving term tP\<p > and the 
rescattering terms in t (NN force contributions only) are peaked at forward angles. Therefore 
the elastic scattering cross section exhibits a characteristic minimum in the angular range 
where the contributions of the exchange and the rescattering terms are of comparable order 
and both are small. This angular range around the minimum could thus be a place where the 
3NF signal, if sufficiently strong, should appear. It would happen at those energies where the 
pure 3NF contribution to the elastic scattering amplitude in that minimum is comparable 
or larger than the contributions of the exchange part and the pure 2N rescattering terms. 

The pure 3NF contribution to the transition operator U results from Eqs. ([I]) and (0) 
when only the 3NF is active: 

jjSNF = pfSNF + yW^ + p) + y(D^ + p^fSNF (4) 

with 

T 3ArF |0 > = F 4 (1) (l + P)\<t> > +F 4 (1) (l + P)G T 3NF \<p > (5) 

We expect that the contribution of U 3NF alone is uniformly distributed over all angles. 

In order to check these expectations we solved Eqs. (0) and ([5]) at the nucleon labora- 
tory energies of 12, 65, 140 and 200 MeV using the modern NN interactions: AV18 jj, CD 
Bonn |3J], Nijm I and Nijm II 0. As the 3NF we took the 27r-exchange Tucson- Melbourne 
(TM) model [[DJ where the strong cut-off parameter A has been adjusted individually to- 



gether with each NN force to the experimental triton binding [|]. In the calculations includ- 
ing 3NF's all partial wave states with total angular momenta in the two-nucleon subsystem 
up to i max = 3 were taken into account. It is the most extensive calculation with 3NF's in 



the continuum which we can perform presently. At the higher energies they are not fully 
converged with respect to j max . The importance of partial waves with higher two-nucleon 
angular momenta are illustrated in fully converged solutions in the case when only 2N forces 
are active. Then we included all states up to j max = 5. Our theoretical results are shown 
in Figs. 1-4 in comparison to data. Our theory does not include the pp Coulomb force. 
Therefore we should compare to nd data. This is only possible at rather low energies, where 
nd data exist and which agree perfectly with NN force predictions only |§. The pd data 
also existing there agree with the nd data except at very forward angles, where Rutherford 
scattering has to show up. That interference with Rutherford scattering can clearly be seen 
in Figs.l and 2 at forward angles, where the data bend towards smaller values. Aside from 
that there is a very good agreement at 12 MeV with theory. This together with the smallness 
of the Coulomb force effects on the elastic scattering cross section in the region of its min- 



imum, as shown by exact calculations under the deuteron breakup threshold |L6]|, supports 
the conjecture that a comparison of nd theory with pd data at even higher energies makes 
sense. The Figs. 1-4 show the expected result, that the pure 3NF contribution is essentially 
uniform in its angular dependence and we see that at 12 MeV it is totally negligible. At 
65 MeV there are also a few nd data [T7J and as shown in Fig. 2 they come close to NN 



force predictions only, whereas the pd data deviate strongly in the minimum. Without 
a rigorous calculation including the pp Coulomb force it has to remain an open question 
whether the deviation between the pd data and the NN force predictions only is due to our 
neglect of Coulomb forces in the theoretical calculations. On the other hand the nd data 
of Fig.3are compatible with pd data in this energy range and indicate only small Coulomb 
force effects corresponding to our conjecture. Apparently precise nd data in the angular 
range of the minima for 65 MeV and higher would be highly desireable. Independent of that 
important issue we can go ahead and display possible 3NF effects in these minima. The 
discrepancy of the theory based on NN forces only to the pd data increases with energy as 
seen in Figs. 2-4. Higher angular momentum states do not cure that discrepancy. They have 
a significant contribution, however, to the cross section at the higher energies at forward 
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angles JT9[, as seen especially in Figs. 3 and 4. As expected the pure 3NF contribution 



remains essentially uniform also at the higher energies. With increasing energy, however, 
this contribution becomes significant in relation to the minimum value of the cross section. 
Being totally negligeable at 12 MeV it overshoots the minimum value by a factor of ~ 6 at 
200 MeV. At 65 MeV the 3NF signal becomes sufficiently large to be seen in the minimum 
region. Indeed as shown in Figs. 2-4, including the 3NF in addition to the 2N interactions in 
the 3N Hamiltonian removes a large part of the discrepancy in the cross section minimum at 
the higher energies. Though the calculations with 3NF are not yet fully converged because of 
presently existing computer limitations we consider that filling of the minima as a "smoking 
gun" for 3NF effects. Very precise data, both in the nd and pd system, would therefore be 
highly valuable. 

We have to expect additional modifications, especially at the highest energies, due to 
relativistic effects, which have not been taken into account in our calculation. First estimates 
just based on kinematical factors indicate indeed a small shift of all angular distribution at 
higher energies toward higher values. 

Finally we want to emphasize that our conclusions do not depend on the particular 
NN interaction used. As shown in Fig. 5 taking different modern NN interactions and the 
corresponding TM 3NF leads to practically the same results. 

Summarizing, we have shown that the minima of the elastic Nd scattering cross sections 
are probably a "smoking gun" for 3NF effects. A large part of the discrepancy between 
modern NN potential predictions and data in this angular range can be removed when the 
TM 3NF properly adjusted to the triton binding is included in the 3N Hamiltonian. In 
order to check more accurately this conclusion precise Nd elastic scattering data at different 
energies in the region of the cross section minima are required. The optimal data would be 
in the nd system to avoid the theoretical uncertainty of pp Coulomb force effects. 
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Figure Captions 



Fig.l The differential Nd cross section at E l ^ b = 12 MeV. The prediction of the CD Bonn 
NN interaction without (short dashed) and with 3NF (solid curve) is compared to pd 



data (circles from [21] and crosses from [22]). The long dashed curve is the pure 3NF 



prediction. All the calculations are truncated at j max = 3. 

Fig. 2 The differential Nd cross section at E l ^ b = 65 MeV. The prediction of the CD Bonn 
NN interaction for j max = 3 (short dashed curve) and j max = 5 (long dashed curve) are 
compared to 64.5 MeV pd data (circles from and nd data (crosses from ]|17[|). 

The CD Bonn calculation including the 3NF for j max = 3 fills the minimum (solid 
curve). The pure 3NF prediction is shown as intermediately long dashed curve. 

Fig. 3 The differential Nd cross section at E l § b = 140 MeV. Curve descritpions as in Fig.2. 



The pd data are: 145.5 MeV (circles) from [24] and 146 MeV (crosses) from p3 
The triangles are 152 MeV nd data from p5 



Fig. 4 The differential Nd cross section at E l £ b = 200 MeV. Curve descritpions as in Fig.2. 
The pd data are: 198 MeV (circles) from @], 200 MeV (crosses) from 0, 181 MeV 



(triangles) from [^3] and 216.5 MeV (x) from ] f2~4[ . 

Fig. 5 The theoretical differential nd cross sections at E l § b = 200 MeV for j max = 3. The 
lower group of curves in the minimum are: AV18 (short dashed), CD Bonn (solid), 
Nijml (intermediately long dashed) and Nijm II (long dashed); the upper group of 
curves in the minimum are the corresponding predictions with TM 3NF included. 
This demonstrates the independence from the specific choice of the NN force. 
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